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T his book has evolved from lectures on engineering mathematics given regu- 
larly over many years to students at all levels in the United States, England, and 
elsewhere. It covers the more advanced aspects of engineering mathematics that 
are common to all first engineering degrees, and it differs from texts with similar 
names by the emphasis it places on certain topics, the systematic development 
of the underlying theory before making applications, and the inclusion of new 
material. Its special features are as follows. 



Prerequisites 

T he opening chapter, which reviews mathematical prerequisites, serves two 
purposes. The first is to refresh ideas from previous courses and to provide 
basic self-contained reference material. The second is to remove from the main 
body of the text certain elementary material that by tradition is usually reviewed 
when first used in the text, thereby allowing the development of more advanced 
ideas to proceed without interruption. 



Worked Examples 

T he numerous worked examples that follow the introduction of each new idea 
serve in the earlier chapters to illustrate applications that require relatively little 
background knowledge. The ability to formulate physical problems in mathemat- 
ical terms is an essential part of all mathematics applications. Although this is not 
a text on mathematical modeling, where more complicated physical applications 
are considered, the essential background is first developed to the point at which 
the physical nature of the problem becomes clear. Some examples, such as the 
ones involving the determination of the forces acting in the struts of a framed 
structure, the damping of vibrations caused by a generator and the vibrational 
modes of clamped membranes, illustrate important mathematical ideas in the 
context of practical applications. Other examples occur without specific applica- 
tions and their purpose is to reinforce new mathematical ideas and techniques as 
they arise. 

A different type of example is the one that seeks to determine the height 
of the tallest flagpole, where the height limitation is due to the phenomenon of 
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buckling. Although the model used does not give an accurate answer, it provides a 
typical example of how a mathematical model is constructed. It also illustrates the 
reasoning used to select a physical solution from a scenario in which other purely 
mathematical solutions are possible. In addition, the example demonstrates how 
the choice of a unique physically meaningful solution from a set of mathematically 
possible ones can sometimes depend on physical considerations that did not enter 
into the formulation of the original problem. 



Exercise Sets 

T he need for engineering students to have a sound understanding of mathe- 
matics is recognized by the systematic development of the underlying theory 
and the provision of many carefully selected fully worked examples, coupled with 
their reinforcement through the provision of large sets of exercises at the ends 
of sections. These sets, to which answers to odd-numbered exercises are listed at 
the end of the book, contain many routine exercises intended to provide practice 
when dealing with the various special cases that can arise, and also more chal- 
lenging exercises, each of which is starred, that extend the subject matter of the 
text in different ways. 

Although many of these exercises can be solved quickly by using standard 
computer algebra packages, the author believes the fundamental mathematical 
ideas involved are only properly understood once a significant number of exer- 
cises have first been solved by hand. Computer algebra can then be used with 
advantage to confirm the results, as is required in various exercise sets. Where 
computer algebra is either required or can be used to advantage, the exercise 
numbers are in blue. A comparison of computer-based solutions with those ob- 
tained by hand not only confirms the correctness of hand calculations, but also 
serves to illustrate how the method of solution often determines its form, and 
that transforming one form of solution to another is sometimes difficult. It is 
the author's belief that only when fundamental ideas are fully understood is it 
safe to make routine use of computer algebra, or to use a numerical package 
to solve more complicated problems where the manipulation involved is pro- 
hibitive, or where a numerical result may be the only form of solution that is 
possible. 



New Material 

T ypical of some of the new material to be found in the book is the matrix 
exponential and its application to the solution of linear systems of ordinary 
differential equations, and the use of the Green's function. The introductory dis- 
cussion of the development of discontinuous solutions of first order quasilinear 
equations, which are essential in the study of supersonic gas flow and in vari- 
ous other physical applications, is also new and is not to be found elsewhere. 
The account of the Laplace transform contains more detail than usual. While 
the Laplace transform is applied to standard engineering problems, including 



control theory, various nonstandard problems are also considered, such as the 
solution of a boundary value problem for the equation that describes the bend- 
ing of a beam and the derivation of the Laplace transform of a function from 
its differential equation. The chapter on vector integral calculus first derives and 
then applies two fundamental vector transport theorems that are not found in 
similar texts, but which are of considerable importance in many branches of 
engineering. 



Series Solutions of Differential Equations 

U nderstanding the derivation of series solutions of ordinary differential equa- 
tions is often difficult for students. This is recognized by the provision of 
detailed examples, followed by carefully chosen sets of exercises. The worked ex- 
amples illustrate all of the special cases that can arise. The chapter then builds 
on this by deriving the most important properties of Legendre polynomials and 
Bessel functions, which are essential when solving partial differential equations 
involving cylindrical and spherical polar coordinates. 



Complex Analysis 

B ecause of its importance in so many different applications, the chapters on 
complex analysis contain more topics than are found in similar texts. In partic- 
ular, the inclusion of an account of the inversion integral for the Laplace transform 
makes it possible to introduce transform methods for the solution of problems 
involving ordinary and partial differential equations for which tables of transform 
pairs are inadequate. To avoid unnecessary complication, and to restrict the mate- 
rial to a reasonable length, some topics are not developed with full mathematical 
rigor, though where this occurs the arguments used will suffice for all practical 
purposes. If required, the account of complex analysis is sufficiently detailed for 
it to serve as a basis for a single subject course. 



Conformal Mapping and Boundary 
Value Problems 

S ufficient information is provided about conformal transformations for them to 
be used to provide geometrical insight into the solution of some fundamen- 
tal two-dimensional boundary value problems for the Laplace equation. Physi- 
cal applications are made to steady-state temperature distributions, electrostatic 
problems, and fluid mechanics. The conformal mapping chapter also provides 
a quite different approach to the solution of certain two-dimensional boundary 
value problems that in the subsequent chapter on partial differential equations 
are solved by the very different method of separation of variables. 



Partial Differential Equations 

A n understanding of partial differential equations is essential in all branches of 
engineering, but accounts in engineering mathematics texts often fall short of 
what is required. This is because of their tendency to focus on the three standard 
types of linear second order partial differential equations, and their solution by 
means of separation of variables, to the virtual exclusion of first order equations 
and the systems from which these fundamental linear second order equations are 
derived. Often very little is said about the types of boundary and initial condi- 
tions that are appropriate for the different types of partial differential equations. 
Mention is seldom if ever made of the important part played by nonlinearity in 
first order equations and the way it influences the properties of their solutions. 
The account given here approaches these matters by starting with first order 
linear and quasilinear equations, where the way initial and boundary conditions 
and nonlinearity influence solutions is easily understood. The discussion of the 
effects of nonlinearity is introduced at a comparatively early stage in the study 
of partial differential equations because of its importance in subjects like fluid 
mechanics and chemical engineering. The account of nonlinearity also includes 
a brief discussion of shock wave solutions that are of fundamental importance in 
both supersonic gas flow and elsewhere. 

Linear and nonlinear wave propagation is examined in some detail because 
of its considerable practical importance; in addition, the way integral transform 
methods can be used to solve linear partial differential equations is described. 
From a rigorous mathematical point of view, the solution of a partial differential 
equation by the method of separation of variables only yields a formal solution, 
which only becomes a rigorous solution once the completeness of any set of 
eigenfunctions that arises has been established. To develop the subject in this 
manner would take the text far beyond the level for which it is intended and 
so the completeness of any set of eigenfunctions that occurs will always be as- 
sumed. This assumption can be fully justified when applying separation of vari- 
ables to the applications considered here and also in virtually all other practical 
cases. 



Technology Projects 

T o encourage the use of technology and computer algebra and to broaden the 
range of problems that can be considered, technology-based projects have 
been added wherever appropriate; in addition, standard sets of exercises of a 
theoretical nature have been included at the ends of sections. These projects are 
not linked to a particular computer algebra package: Some projects illustrating 
standard results are intended to make use of simple computer skills while others 
provide insight into more advanced and physically important theoretical ques- 
tions. Typical of the projects designed to introduce new ideas are those at the 
end of the chapter on partial differential equations, which offer a brief introduc- 
tion to the special nonlinear wave solutions called solitons. 



Numerical Mathematics 



A lthough an understanding of basic numerical mathematics is essential for all 
engineering students, in a book such as this it is impossible to provide a sys- 
tematic account of this important discipline. The aim of this chapter is to provide 
a general idea of how to approach and deal with some of the most important 
and frequently encountered numerical operations, using only basic numerical 
techniques, and thereafter to encourage the use of standard numerical packages. 
The routines available in numerical packages are sophisticated, highly optimized 
and efficient, but the general ideas that are involved are easily understood once 
the material in the chapter has been assimilated. The accounts that are given 
here purposely avoid going into great detail as this can be found in the quoted 
references. However, the chapter does indicate when it is best to use certain types 
of routine and those circumstances where routines might be inappropriate. 

The details of references to literature contained in square brackets at the ends 
of sections are listed at the back of the book with suggestions for additional read- 
ing. An instructor's Solutions Manual that gives outline solutions for the techno- 
logy projects is also available. 
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PART ONE 



REVIEW MATERIAL 




CHAPTER 



S' 



Review of Prerequisites 



E very account of advanced engineering mathematics must rely on earlier mathematics 
courses to provide the necessary background. The essentials are a first course in calculus 
and some knowledge of elementary algebraic concepts and techniques. The purpose of 
the present chapter is to review the most important of these ideas that have already been 
encountered, and to provide for convenient reference results and techniques that can be 
consulted later, thereby avoiding the need to interrupt the development of subsequent 
chapters by the inclusion of review material prior to its use. 

Some basic mathematical conventions are reviewed in Section 1 .1, together with the 
method of proof by mathematical induction that will be required in later chapters. The 
essential algebraic operations involving complex numbers are summarized in Section 1 .2, 
the complex plane is introduced in Section 1.3, the modulus and argument representa- 
tion of complex numbers is reviewed in Section 1.4, and roots of complex numbers are 
considered in Section 1 .5. Some of this material is required throughout the book, though 
its main use will be in Part 5 when developing the theory of analytic functions. 

The use of partial fractions is reviewed in Section 1 .6 because of the part they play 
in Chapter 7 in developing the Laplace transform. As the most basic properties of deter- 
minants are often required, the expansion of determinants is summarized in Section 1.7, 
though a somewhat fuller account of determinants is to be found later in Section 3.3 of 
Chapter 3. 

The related concepts of limit, continuity, and differentiability of functions of one or 
more independent variables are fundamental to the calculus, and to the use that will 
be made of them throughout the book, so these ideas are reviewed in Sections 1 .8 and 
1.9. Tangent line and tangent plane approximations are illustrated in Section 1.10, and 
improper integrals that play an essential role in the Laplace and Fourier transforms, and 
also in complex analysis, are discussed in Section 1.11. 

The importance of Taylor series expansions of functions involving one or more in- 
dependent variables is recognized by their inclusion in Section 1.12. A brief mention is 
also made of the two most frequently used tests for the convergence of series, and of the 
differentiation and integration of power series that is used in Chapter 8 when considering 
series solutions of linear ordinary differential equations. These topics are considered again 
in Part 5 when the theory of analytic functions is developed. 

The solution of many problems involving partial differential equations can be simplified 
by a convenient choice of coordinate system, so Section 1 .1 3 reviews the theorem for the 
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change of variable in partial differentiation, and describes the cylindrical polar and spherical 
polar coordinate systems that are the two that occur most frequently in practical problems. 

Because of its fundamental importance, the implicit function theorem is stated without 
proof in Section 1 .14, though it is not usually mentioned in first calculus courses. 



1.1 Real Numbers, Mathematical induction, 

m and Mathematical Conventions 

N umbers are fundamental to all mathematics, and real numbers are a subset 
of complex numbers. A real number can be classified as being an integer, a 
rational number, or an irrational number. From the set of positive and negative 
integers, and zero, the set of positive integers 1, 2, 3, . . . is called the set of natural 
numbers. The rational numbers are those that can be expressed in the form m/n, 
where m and n are integers with n / 0. Irrational numbers such as n, ~J2, and sin 2 
are numbers that cannot be expressed in rational form, so, for example, for no 
integers m and n is it true that \fl is equal to m/n. Practical calculations can only 
be performed using rational numbers, so all irrational numbers that arise must be 
approximated arbitrarily closely by rational numbers. 

Collectively, the sets of integers and rational and irrational numbers form what 
is called the set of all real numbers, and this set is denoted by R. When it is necessary 
to indicate that an arbitrary number a is a real number a shorthand notation is 
adopted involving the symbol e, and we will write a e R. The symbol e is to be read 
“belongs to” or, more formally, as “is an element of the set.” If a is not a member 
of set R, the symbol e is negated by writing </. and we will write a </ R where, of 
course, the symbol is to be read as “does not belong to,” or “is not an element 
of the set.” As real numbers can be identified in a unique manner with points on a 
line, the set of all real numbers R is often called the real line. The set of all complex 
numbers C to which R belongs will be introduced later. 

One of the most important properties of real numbers that distinguishes them 
from other complex numbers is that they can be arranged in numerical order . This 
fundamental property is expressed by saying that the real numbers possess the order 
property. This simply means that if x, y e R, with x / y, then 

either x < y or x > y, 

where the symbol < is to be read “is less than” and the symbol > is to be read 
“is greater than.” When the foregoing results are expressed differently, though 
equivalently, if x, y e R, with x ^ y, then 

either x — y < 0 or x — y > 0. 

It is the order property that enables the graph of a real function / of a real 
variable x to be constructed. This follows because once length scales have been 
chosen for the axes together with a common origin, a real number can be made 
to correspond to a unique point on an axis. The graph of / follows by plotting all 
possible points (x, f(x)) in the plane, with x measured along one axis and f(x) 
along the other axis. 

The absolute value |jc| of a real number x is defined by the formula 



absolute value 



x 



x if x > 0 
—x if x < 0. 
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This form of definition is in reality a concise way of expressing two separate state- 
ments. One statement is obtained by reading |x with the top condition on the right 
and the other by reading it with the bottom condition on the right. The absolute 
value of a real number provides a measure of its magnitude without regard to its 
sign so, for example, |3| = 3, | — 7.41 1 = 7.41, and |0| = 0. 

Sometimes the form of a general mathematical result that only depends on an 
arbitrary natural number n can be found by experiment or by conjecture, and then 
the problem that remains is how to prove that the result is either true or false for 
all n. A typical example is the proposition that the product 

(1 - 1/4)(1 - 1/9)(1 - 1/16) ... [1 - l/(« + l) 2 ] 

= (n + 2)/(2n + 2), for n = 1,2 

This assertion is easily checked for any specific positive integer n, but this does not 
amount to a proof that the result is true for all natural numbers. 

A powerful method by which such propositions can often be shown to be either 
true or false involves using a form of argument called mathematical induction. This 
type of proof depends for its success on the order property of numbers and the fact 
that if n is a natural number, then so also is n + 1. The steps involved in an inductive 
proof can be summarized as follows. 

Proof by Mathematical Induction 

Let P(n) be a proposition depending on a positive integer n. 

STEP 1 Show, if possible, that P(n) is true for some positive integer no. 

STEP 2 Show, if possible, that if P(n) is true for an arbitrary integer n = k> no, 

then the proposition P(k+ 1) follows from proposition P(k). 

STEP 3 if Step 2 is true, the fact that P(n o) is true implies that P(?io + 1) is true, 
and then that P{uq + 2) is true, and hence that P(n) is true for all n > no. 
STEP 4 If no number n = no can be found for which Step 1 is true, or if in Step 2 
it can be shown that P{k) does not imply P(k +1), the proposition P(n) 
is false. 

The example that follows is typical of the situation where an inductive proof 
is used. It arises when determining the nth term in the Maclaurin series for sin ax 
that involves finding the nth derivative of sin ax. A result such as this may be found 
intuitively by inspection of the first few derivatives, though this does not amount to 
a formal proof that the result is true for all natural numbers n. 

Prove by mathematical induction that 

d n /dx" [sin ax] = a n sin(ax + nit IT), for n — 1,2, 

Solution The proposition P(n) is that 

d n /dx n [s\nax] = a n sin(ax + mr/2), for n — 1,2, 

STEP I Differentiation gives 



d/dx\sinax ] = a cosax, 
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but setting n— 1 in P(n) leads to the result 

d /dx[smax\ — a sin(nx + n/2) = a cos ax, 



showing that proposition P{n) is true for n = 1 (so in this case n 0 = 1). 
STEP 2 Assuming P(k) to be true for k > 1, differentiation gives 



d/dx{d k /dx k [smax]} — d/dx[a k sin(ax + for/2)], 



so 



d k+1 /dx k+1 [sinax] = a k+l cos {ax + kjt/2). 



However, replacing k by k + 1 in P(k) gives 

d k+1 / dx k+1 [sin ax] = a k+l sin[n.r + (k+ l)?r/2] 

= a k+1 sin[(n.r + kjt/2) + 7r/2] 

= a k+l cos {ax + for/2), 

showing, as required, that proposition P(k ) implies proposition P(k+ 1), so Step 2 
is true. 

STEP 3 As P(n) is true for n = 1, and P(k ) implies P(k+ 1), it follows that the 
result is true for n — 1 , 2 ,... and the proof is complete. 

The binomial theorem finds applications throughout mathematics at all levels, 
so we quote it first when the exponent n is a positive integer, and then in its more 
general form when the exponent a involved is any real number. 

Binomial theorem when n is a positive integer 

If a, b are real numbers and n is a positive integer, then 



( fl + b) n = a n + na n ~ l b + 




n(n — 1 ){n — 2) 
+ 3! 



a 



,»- 3 1,3 



b 3 + --- + b'\ 



binomial coefficient 



or more concisely in terms of the binomial coefficient 




we have 




r = 0 



where m\ is the factorial function defined as m\ = 1 • 2 ■ 3 • • • m with m > 0 an 
integer, and 0! is defined as 0! = 1. It follows at once that 
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multi-statement 



The binomial theorem involving the expression ( a + b) a , where a and b are real 
numbers with \b/a\ < 1 and a is an arbitrary real number takes the following form. 



General form of the binomial theorem when a is an arbitrary 
real number 

If a and b are real numbers such that \b/a\ < 1 and a is an arbitrary real 
number, then 



(a + bf = a 01 ( 1 + - 

V a 




a(u — l)/b 
2! Vfl 



2 



a(a — l)(oi — 2) 

+ 3! 




The series on the right only terminates after a finite number of terms if a is a 
positive integer, in which case the result reduces to the one just given. If a is 
a negative integer, or a nonintegral real number, the expression on the right 
becomes an infinite series that diverges if \b/a\ > 1. 



Expand (3 + x) 1/2 by the binomial theorem, stating for what values of .r the series 
converges. 

Solution Setting b/a = \x in the general form of the binomial theorem gives 

(3 - *)-** = r*(i - + 2 V - + ' ' O' 

The series only converges if \l x \ < !> and so it is convergent provided |x| < 3. ■ 



Some standard mathematical conventions 

Use of combinations of the ± and ^ signs 

The occurrence of two or more of the symbols ± and =p in an expression is to be 
taken to imply two separate results, the first obtained by taking the upper signs and 
the second by taking the lower signs. Thus, the expression a ± b sin 6 =p c cos 6 is an 
abbreviation for the two separate expressions 

a + bsin9 — ccost? and a — bsin9 + ccost?. 

Multi-statements 

When a function is defined sectionally on n different intervals of the real line, instead 
of formulating n separate definitions these are usually simplified by being combined 
into what can be considered to be a single multi-statement. The following example 
is typical of a multi-statement: 

sinx, x < it 
0, n < x < 3 tt /2 
— 1, x > 3tt/2. 



/(*) = 
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It is, in fact, three statements. The first is obtained by reading f(x ) in conjunction 
with the top line on the right, the second by reading it in conjunction with the second 
line on the right, and the third by reading it in conjunction with the third line on 
the right. An example of a multi-statement has already been encountered in the 
definition of the absolute value |jc| of a number x. Frequent use of multi-statements 
will be made in Chapter 9 on Fourier series, and elsewhere. 

Polynomials 

A polynomial is an expression of the form P(x) = a^x" + ayx n ~ 1 -\- ■ ■ ■ + a n ~ ix + a n ■ 
The integer n is called the degree of the polynomial, and the numbers a t are called 
its coefficients. The fundamental theorem of algebra that is proved in Chapter 14 
asserts that P(x) = 0 has n roots that may be either real or complex, though some 
of them may be repeated, (no ^ 0 is assumed.) 

Notation for ordinary and partial derivatives 

If f(x) is an n times differentiable function then f (n \x) will, on occasion, be used 
to signify d n f/dx n , so that 



/ (n) « 



d^l 

dx n ' 



suffix notation for 
partial derivatives 



If f(x, y) is a suitably differentiable function of x and y, a concise notation used to 
signify partial differentiation involves using suffixes, so that 



3 / 



3 ( 3 / 



fx = — , fyx = (fy)x = 7- h— 



dx 



3 2 f 



dx \d y ) dydx 



f -tl 

hy ~ dy 2 ’ 



with similar results when / is a function of more than two independent variables. 
Inverse trigonometric functions 

The periodicity of the real variable trigonometric sine, cosine, and tangent functions 
means that the corresponding general inverse trigonometric functions are many val- 
ued. So, for example, if y = sin x and we ask for what values of x is y = 1 /s/2, we 
find this is true for x = jt/4 ± 2 nix and x = 3jz/4 ± 2mt for n = 0, 1, 2 To over- 

come this ambiguity, we introduce the single valued inverses, denoted respectively 
by x — Arcsin y, x = Arccos y, and x = Arctan y by restricting the domain and 
range of the sine, cosine, and tangent functions to one where they are either strictly 
increasing or strictly decreasing functions, because then one value of x corresponds 
to one value of y and, conversely, one value of y corresponds to one value of x. 

In the case of the function y = sin x, by restricting the argument x to the inter- 
val — 7t/2 < x < tt/2 the function becomes a strictly increasing function of x. The 
corresponding single valued inverse function is denoted by x = Arcsin y, where y is 
a number in the domain of definition [—1, 1] of the Arcsine function and x is a num- 
ber in its range [—tt/2, tt/2]. Similarly, when considering the function y = cos x, the 
argument is restricted to 0 < x < n to make cos x a strictly decreasing function of x. 
The corresponding single valued inverse function is denoted by x = Arccos y, where 
y is a number in the domain of definition [—1,1] of the Arccosine function and x is 
a number in its range [0, n]. Finally, in the case of the function y = tan x, restricting 
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the argument to the interval —tt/2 < x < tt / 2 makes the tangent function a strictly 
increasing function of x. The corresponding single valued inverse function is de- 
noted by x = Arctan y where y is a number in the domain of definition (— oo, oo) 
of the Arctangent function and a: is a number in its range {—tt/2, tt/2). 

As the inverse trigonometric functions are important in their own right, the 
variables x and y in the preceding definitions are interchanged to allow considera- 
tion of the inverse functions y = Arcsin x , y = Arccos x, and y = Arctan x, so that 
now x is the independent variable and y is the dependent variable. 

With this interchange of variables the expression y = arcsin x will be used to 
refer to any single valued inverse function with the same domain of definition as 
Arcsin x, but with a different range. Similar definitions apply to the functions y = 
arccos x and y = arctan x. 

Double summations 

An expression involving a double summation like 

OO OO 

EE a mn sin mx sin ny , 

m= 1 n= 1 

means sum the terms a mn sin mx sin ny over all possible values of m and n, so that 

OO OO 

EE a mn sin mx sin ny — an sin x sin y + a u sin x sin 2y 

m= 1 n = 1 

+ « 2 i sin 2 x sin y + a 2 2 sin 2 x sin 2 y H . 

A more concise notation also in use involves writing the double summation as 

OO 

E a mn sin rax sin nx. 

m=l,n=l 



signum function 



The signum function 



The signum function, usually written sign(x), and sometimes sgn(x), is defined as 



sign(x) = 




if x > 0 
if x < 0. 



We have, for example, sign(cosx) = 1 for 0 < x < Tt/2, and sign(cosx) = — 1 for 
7r/2 < x < tt or, equivalently, 



sign(cos x) 



1, 0 < x < ^ 7 r 

— 1, < X < Tt. 



infinite product 



Products 

Let {iik\\ \-i be a sequence of numbers or functions u\, w 2 , . . . ; then the product of 
the n members of this sequence is denoted by n/=i M *> so that 

n 

~ | life = U l //? ■ ■ ■ u n . 
k= 1 



iim n ii k = r 

n — >oo 1 1 



U k 



k=l k = 1 



When the sequence is infinite, 
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is called an infinite product involving the sequence {iik}. Typical examples of infinite 
products are 





- and 
2 




sin x 
x 



More background information and examples can be found in the appropriate sec- 
tions in any of references [1.1], [1.2], and [1.5]. 



Logarithmic functions 

The notation In x is used to denote the natural logarithm of a real number x. that 
is, the logarithm of x to the base e, and in some books this is written log e x. In this 
book logarithms to the base 10 are not used, and when working with functions of a 
complex variable the notation Log z, with z = re tB means Log z = lnr + id. 



EXERCISES 1.1 



1. Prove that if a > 0, b > 0, then a/Vb + b/ fa > 

fa + Vb. 

Prove Exercises 2 through 6 by mathematical induction. 

2. YZkZo( a + kd) = (n/2)[2a + ( n - l)rf] 

(sum of an arithmetic series). 

3- E^Jr* = (l-r")/(l-r) (r * 1) 

(sum of a geometric series). 

4. YTk=\ & — (1/6 )n(n + l)(2n + 1) (sum of squares). 

5. d n /dx n [cosax] = a n cos(ax + mc/2), with n a natural 
number. 

6. d n /dx n [ ln(l + x)\ = (— 1 )" +1 (« — 1)!/(1 + x) n , with n a 
natural number. 



7. Use the binomial theorem to expand (3 + 2x) 4 . 

8. Use the binomial theorem and multiplication to expand 
(1 —x 2 )(2 + 3x) 3 . 

In Exercises 9 through 12 find the first four terms of the 
binomial expansion of the function and state conditions for 
the convergence of the series. 

9. (3 + 2x)~ 2 . 

10. (2 — x 2 ) 1/3 . 

11. (4 + 2 x 2 )- 1 / 2 . 

12. (1 - 3x 2 ) 3/4 . 




Complex Numbers 

Mathematical operations can lead to numbers that do not belong to the real number 
system R introduced in Section 1.1. In the simplest case this occurs when finding 
the roots of the quadratic equation 



ax 1 + bx + c = 0 with a, b, c e R, a ^ 0 



by means of the quadratic formula 

—b ± Vb 2 — 4 ac 



The discriminant of the equation is b 2 — 4 ac, and if b 2 — Aac < 0 the formula 
involves the square root of a negative real number; so, if the formula is to have 
meaning, numbers must be allowed that lie outside the real number system. 

The inadequacy of the real number system when considering different math- 
ematical operations can be illustrated in other ways by asking, for example, how 
to find the three roots that are expected of a third degree algebraic equation as 
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real and imaginary 
part notation 



simple as x 3 — 1 = 0, where only the real root 1 can be found using y = x 3 — 1, or 
by seeking to give meaning to ln(— 1), both of which questions will arise later. 

Difficulties such as these can all be overcome if the real number system is 
extended by introducing the imaginary unit i defined as 

i 2 = -1, 

so expressions like -y/(— k 2 ) where k a positive real number may be written 
v /( — 1 ) v /(/c 2 ) = ±ik. Notice that as the real number k only scales the imaginary 
unit i, it is immaterial whether the result is written as ik or as ki. 

The extension to the real number system that is required to resolve problems 
of the type just illustrated involves the introduction of complex numbers, denoted 
collectively by C, in which the general complex number, usually denoted by z„ has 
the form 



z — a + ip, with a, ft real numbers. 

The real number a is called the real part of the complex number z, and the real 
number p is called its imaginary part. When these need to be identified separately, 
we write 



Re{z} = a and lm{z} = p, 

so if z = 3 — li, Re{z) = 3 and Im{z} = —7. 

If Im{z} — P = 0 the complex number z reduces to a real number, and if 
Re{z} = a = 0 it becomes a purely imaginary number, so, for example, z = 5 i is 
a purely imaginary number. When a complex number z is considered as a variable 
it is usual to write it as 

Z = x + iy, 

where x and y are now real variables. If it is necessary to indicate that z is a general 
complex number we write z e C. 

When solving the quadratic equation az 2 + bz + c = 0 with a, b, and c real 
numbers and a discriminant b 2 — 4a c < 0, by setting 4 ac — b 2 = k 2 in the quadratic 
formula, with k> 0, the two roots zi and zi are given by the complex numbers 

Zi = —( b/2a ) + i(k/2a) and zi = —(b/2a) — i(k/2a). 

Algebraic rules for complex numbers 

Let the complex numbers zi and z .2 be defined as 

zi=a+ib and zi = c + id, 

with a,b,c, and d arbitrary real numbers. Then the following rules govern the 
arithmetic manipulation of complex numbers. 



Equality of complex numbers 

The complex numbers z\ and z 2 are equal, written z\ = z 2 if, an d only if, 
Rejzi) = Re{z 2 } and Im{zi) = Im{z 2 ). So a + ib = c + id if, and only if, 

a = c and b = d. 
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EXAMPLE 1.4 



EXAMPLE 1.5 



(a) 3 — 9i = 3 + bi if, and only if, b = —9. 

(b) If u = — 2 + 5i, v = 3 + 5 i, w = a + 5 i, then 
u = w if, and only if, a = — 2 but u ^ v, and 

v = w if, and only if, a = 3. ■ 

Zero complex number 

The zero complex number, also called the null complex number, is the number 
0 + 0 i that, for simplicity, is usually written as an ordinary zero 0. 



If a + ib = 0, then a = 0 and b = 0. 

Addition and subtraction of complex numbers 

The addition (sum) and subtraction (difference) of the complex numbers zi 
and zz is defined as 

zi + Z 2 = Refzi) + Re{z 2 } + i[Im{zi} + Im{z 2 }] 

and 

Z] -z,2 = Re{zi} - Re{z 2 ) + /' [Im{zi } - Im{z 2 }]. 



So, if zi = a + ib and z 2 = c + id, then 



and 



Zi + Z2 = (a + ib) + (c + id) 
= (a + c) + i(b + d), 

Z\- Z2 = (a + ib) - (c + id) 
= (a — c) + i{b — d). 



If zi = 3 + 7/ and z 2 = 3 + 2i, then the sum 

Zi + Z 2 = (3 + 3) + (7 + 2 )i = 6 + 9 i, 

and the difference 

zi - Z2 = (3 - 3) + (7 - 2 )i = 5 i. 



Multiplication of complex numbers 

The multiplication (product) of the two complex numbers z\ = a + ib and 
Z2 = c + id is defined by the rule 

Z\Z2 — (a + ib)(c + id) = ( ac — bd) + i(ad + be). 



An immediate consequence of this definition is that if A: is a real number, then 
kz\ = k(a + ib) = ka + ikb. This operation involving multiplication of a complex 
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number by a real number is called scaling a complex number. Thus, if zi = 3 + li 
and zi = 3 + 2 i, then 2zi — 3z2 = (6 + 14/) — (9 + 6 i) = — 3 + 8 i. 

In particular, if z = a + ib, then — z — (— l)z — — a — ib. This is as would be 
expected, because it leads to the result z — z = 0. 

In practice, instead of using this formal definition of multiplication, it is more 
convenient to perform multiplication of complex numbers by multiplying the brack- 
eted quantities in the usual algebraic manner, replacing every product i 2 by —1, and 
then combining separately the real and imaginary terms to arrive at the required 
product. 



(a) 5i(-4 + 3i) = —15 — 201 

(b) (3 - 2z')( — 1 + 4i)(l + 0 = (-3 + 12 i + 2 i - 8/ 2 )(l + i) 

= [(— 3 + 8) + (12 + 2)i](l + 0 = (5 + 140(1 + 0 
= 5 + 14/ + 5 i + 14/ 2 = (5 - 14) + (5 + 14)/ = -9 + 19/. 



Complex conjugate 

If z = a + ib, then the complex conjugate of z, usually denoted by z and read 
“z bar,” is defined as z. = a — ib. It follows directly that 

(z) = z and zz = a 2 + b 2 . 



In words, the complex conjugate operation has the property that taking the complex 
conjugate of a complex conjugate returns the original complex number, whereas 
the product of a complex number and its complex conjugate always yields a real 
number. 

If z— a + ib, then adding and subtracting z and z gives the useful results 
z+z = 2Re{z} = 2a and z - z = 2/ Im{z} = 2ib. 



These can be written in the equivalent form 



Re{z} = a = i(z+ z) and Im{z} = b — ^(z - z). 
2 2 1 



Quotient (division) of complex numbers 

Let zi = a + ib and zi = c + id. Then the quotient Z\/Z 2 is defined as 



Zi 

Z2 



(ac + bd) + i{bc — ad) 
c 2 + d 2 



, Z2 ^ 0 . 





